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Abstract 

Oblique propagation of magnetohydrodynamic waves in warm plasmas is described 
by a modified vector derivative nonlinear Schrodinger equation, if charge separation 
in Poisson's equation and the displacement current in Ampere's law are properly 
taken into account. This modified equation cannot be reduced to the standard 
derivative nonlinear Schrodinger equation and hence its possible integrability and 
related properties need to be established afresh. Indeed, the new equation is shown 
to be integrable by the existence of a bi-Hamiltonian structure, which yields the 
recursion operator needed to generate an infinite sequence of conserved densities. 
Some of these have been found explicitly by symbolic computations based on the 
symmetry properties of the new equation. Since the new equation includes as a 
special case the derivative nonlinear Schrodinger equation, the recursion operator 
for the latter one is now readily available. 



1 Introduction 

The derivative nonlinear Schrodinger equation (DNLS) was first given by Rogister (1971) 
for the nonhnear evolution of parallel Alfven waves in plasmas, and later encountered 
in many different contexts by other authors, emerging as one of the canonical nonlinear 
equations in physics. The DNLS could also account for slightly oblique propagation of 
Alfven waves (see e.g. Hada, Kennel and Buti 1989), albeit at the price of neglecting 
two effects which might be important in strongly magnetized astrophysical plasmas. One 
is the deviation from charge neutrality between the different plasma species, the other 
is the influence of the displacement current in Ampere's law. Retaining these effects 
results in a nonlinear vector evolution equation which differs from the standard vector 
form of the DNLS by an extra linear term, and therefore was called the modifled vector 
derivative nonlinear Schrodinger equation (MVDNLS) (Deconinck, Meuris and Verheest 
1993a, 1993b). 

Recently, there has been a renewed interest in the use of the DNLS for certain astro- 
physical plasmas (Spangler 1992, Spangler and Plapp 1992), assuming that the case of 
slightly oblique propagation could easily be reduced to that of parallel propagation. To 
do so, one modifles the dependent variable, in this case the perpendicular magnetic fleld, 
by including the static part as well. Since the MVDNLS cannot be transformed into the 
DNLS itself, such an easy transition from parallel to oblique propagation is not possible. 
The reverse is true, of course, the MVDNLS includes the DNLS as a special case, when 
we go from oblique to parallel propagation, which amounts to dropping the bothersome 
extra term. 

That term, which distinguishes the MVDNLS from the DNLS, has implications for 
the discussion of integrability and the possibility of deriving solitary wave solutions for the 
MVDNLS. The DNLS is well known to be completely integrable (Kaup and Newell 1978), 
whereas for the MVDNLS we could only get certain indications about its integrability 
(Deconinck, Meuris and Verheest 1993b). The applicability of the prolongation method 
(Kaup 1980), adapted to a vector nonlinear equation, and the existence of some invariants 
(Deconinck, Meuris and Verheest 1993b) were indicative of complete integrability, without 
giving a watertight proof. 

In the present paper we show that the MVDNLS possesses a bi-Hamiltonian struc- 
ture, and hence through the resulting recursion operator an inflnite sequence of conserved 
densities. Interestingly enough, to the best of our knowledge there seems to be no proof 
in the literature that the DNLS itself has a bi-Hamiltonian structure, although this ex- 
pected property is mentioned sometimes without further details nor references (Oevel and 



Fuchssteiner 1992). The existence of an infinite sequence of conserved densities is proved 
by Kaup and Newell (1978) for the DNLS, without showing the bi-Hamiltonian charac- 
ter. Our constructions, formulas and conclusions for the MVDNLS immediately apply to 
the DNLS. Hence, without extra work, the explicit forms of both Hamiltonians and the 
recursion operator for the DNLS are now available. 

In §2 we recall the form of the MVDNLS and list in §3 some of the conserved densities, 
which were found in an ad hoc fashion with the help of a symbolic program and by looking 
at the symmetry properties of the equation. The knowledge of these conserved densities 
turned out to be beneficial for an easy construction of the appropriate Hamiltonians in 
§4. In §5 the bi-Hamiltonian structure and the recursion operator are derived and with 
it we established the existence of an infinite sequence of conserved densities, needed to 
guarantee complete integrability. §6 is devoted to a short discussion of the implications 
for the usual DNLS, and in §7 we draw some conclusions. 



2 MVDNLS 

The MVDNLS is, after the necessary scaling and Galilean transforms to cast it in its 
simplest dimensionless form, given by 

where the parameter a characterizes the extra term which distinguishes the MVDNLS 
from the DNLS. In (1) B^ stands for the perpendicular magnetic field, which includes both 
the wave contributions and the static perpendicular field due to the oblique propagation 
with respect to the total external magnetic field. The direction of wave propagation is 
along the x-axis. If the third term is absent, we can project (1) onto axes perpendicular to 
the direction of wave propagation, introduce a new complex variable from the components 
ofB^, 

0± = 5,±2S„ (2) 

and combine the projections to obtain the DNLS in standard scalar form. 



dt dx dx"^ 

The ± signs in (2) and (3) are correlated. The DNLS can account for oblique propa- 
gation, provided a is zero, i.e. imposing charge neutrality to all orders and neglecting 
the displacement current in Ampere's law (Deconinck, Meuris and Verheest 1993a). The 



bothersome third term in (1) also disappears if B_|_o is zero, in the case of strictly parallel 
propagation. 

At this stage, it is worth recalling that the DNLS has constant-amplitude solutions 
of the form 

B,, = acos(kx — ut), 

(4) 

Bz = ± asm{kx — ut), 

with a an arbitrary constant. As can easily be checked, there are no constant amplitude 
solutions to (1) besides the trivial case with B_l = B_lo, if B_lo 7^ 0. This means that any 
sort of separation into left and right circularly polarized waves as for the DNLS is doomed 
to fail. Of course, for oblique propagation with a = 0, the circular polarization given by 
the DNLS is only apparent, since B_|_ includes the static part B_|_o, and this shift leads 
in reality to elliptical polarization for the perpendicular wave field (Spangler and Plapp 
1992). In addition, the MVDNLS has a class of stationary solitary wave solutions which 
the DNLS does not have, the subalfvenic modes, which have totally different properties 
compared to the known stationary solutions of the DNLS (Deconinck, Meuris and Verheest 
1993b). 

We know that the DNLS is completely integrable in the sense that it is possesses an 
infinite series of conserved densities, which can be constructed explicitly. As the MVDNLS 
reduces for a — ^ to the vector form of the DNLS, there is hope to prove integrability 
for the MVDNLS. There are, however, two major complications with the MVDNLS: it 
has a vector character which contrary to the DNLS cannot be transformed away, and the 
boundary values at infinity are not zero (at least not in the physical model for which the 
nonlinear evolution equation was derived!). 

One encounters in the literature quite a variety of methods to investigate symmetries, 
to construct conservation laws, or to establish integrability of nonlinear equations via 
direct or inverse methods. While in principle these methods could be used, one rarely 
sees worked examples involving vector equations. Furthermore, the nonzero boundary 
conditions for the MVDNLS are an additional hurdle. 



3 Integrability and conserved densities 

The first step in any treatment is to try to prove integrability, or at least collect sufficiently 
compelling evidence. For instance, one could check if the equation passes the Painleve test. 



However, here again the nonscalar character of the MVDNLS prevents straightforward 
apphcation of this otherwise so useful test (Fordy 1990). A different way to ascertain 
integrabihty, at least as convincing, is indicated by Kaup (1980) and Fordy (1990). It 
involves the application of the prolongation method due to Estabrook and Wahlquist. 
This method was adapted successfully to the nonlinear equation at hand (Deconinck, 
Meuris and Verheest 1993b). 

Yet another important tool to determine integrabihty of a nonlinear PDE is finding 
a sufficiently large number of conservation laws of the form 

where J^^ p dx is the conserved quantity with density p and associated flux J. We have 
assumed functions u and v which are fast decreasing at inflnity to symmetric but in the 
case of u non-zero boundary conditions. As is well known from many other examples, 
the flrst conservation law comes from rewriting the equations themselves in the form of a 
conservation law. For the MVDNLS this yields after projection 

ut + (u{u^ + t;2) +(3u- v^) = 0, 

where u and v denote the components of B_|_ parallel and perpendicular to B_|_o, and 
f3 = <yBj_Q. As usual, subscripts refer to partial derivatives with respect to t and x. 
Whereas (6) amounts to a vector conservation law, the other ones we have derived through 
constructive procedures are all scalar ones, involving powers of {u^ + f^), as presented in 
(9)-(13). 

To see how to proceed, we follow ideas exposed in more detail in Verheest and Hereman 
(1994), and briefly discuss the scaling or symmetry properties of the equations. These 
can be used to obtain information about polynomial conserved densities and the building 
blocks they are made off. The scaling of (6) is such that 

ur^v, ^ ~ ^-^' M ^v ~/y~--. (7) 

ot ax^ ox 

We may restrict ourselves to building blocks which belong to the same class under the 
mentioned scaling, since for any mixed conserved quantity which one could derive, the 
freedom implied in the scaling would split that quantity in several conserved quantities, 
each with building blocks of the same scaling. Thus there is a straightforward and logical 
way to construct invariant quantities (Miura, Gardner and Kruskal (1968)). Moreover 
and without loss of generality, we may remove any density (or part thereof) that is a 



total x-derivative, for these are trivially conserved. In addition, for /3 = (6) is invariant 
under the substitution 

M — > t;, V — > — -u. (8) 

In every conserved density the part without factors j3 will have to obey this additional 
symmetry. 

At the quadratic level in u and v, the only possibility is u^ + v"^ or (9) given below, 
due to the rule in (8). Cubic terms in u and v are not possible at all. 

Starting then with a candidate density containing the building block (-u^ + f ^)^, one 
has four factors u or v, and one could add a combination of the form u^v — uv^. Keeping 
the scaling (7) in mind, quartic terms are also equivalent to two factors u and/or v and 
one derivation, or to two factors u and/or v and one factor (3. Two factors u and/or v with 
one derivation can only lead to a non-trivial building block of the structure uvx — vux, if 
we keep all the preceding remarks in mind. With one factor 13 we could in principle have 
a linear combination of u^, uv and v"^. This exhausts the building blocks at this order 
and it is then for the MATHEMATICA program we wrote to determine the necessary 
coefficients. This leads to (10). Obviously, we can go on like this for higher orders. 

Labelling the conserved densities p„, with n the corresponding power of {v^ + f ^), we 
obtained 

Pi = u^ + v\ (9) 

X 

P2 = -{U^ + V^y - UVx + u^v + (3u^ , (10) 

P3 = \{u' + vr+^-{ul + vl)-uhx+v\x + (^{u'-v'), (11) 

12 4 

p4 = -:{U^ + V^y - -{UxVxx-UxxVx) + -{uux + vvxf 
4 5 5 

+ -{u^ + v^)iul + vl) - (u^ + v^iuvx - UxV) (12) 

5 

+ ^{2ul - Auhx + 2v^ + 3nS2 - t;6) + ^u\ 
5 5 

P5 = T7:{u'^ + v'^f + T^i'^lx + ^D-o(^^ + v'^){uxVxx-UxxVx) + 5(m^ + v'^){uux + VVxf 

+ -r{u^ + v^f{ul + vlf - —{u" + v''f{uvx - Uxv) 

4 16 (13) 

+ f (5u^ + 10m%2 - 10m\^ - hv^ + 20m\^ - l2u^Vx + muv^Vx - 20v\^) 
8 

4 ^ ' 



Note that we have not yet included the renormahzation constants needed to ensure the 
boundedness of the conserved quantities obtained from the above hsted densities. We will 
come back to this point in the following paragraph. 



4 Hamiltonian structure 



Let us start by pointing out some of the principal ingredients of the Hamiltonian structure 
of evolution equations, adapted where necessary to vector quantities. 

The system (6) is said to possess a Hamiltonian structure (Olver 1980), if there exists 
a so called Hamiltonian operator O (Fokas 1987) (sometimes called implectic operator, see 
e.g. Fuchssteiner and Fokas 1981) and a (2-component) gradient vector function ^^{u^v) 
such that (6) can be written in the form: 




e-7^. 



(14) 



The operator is a Hamiltonian operator if it is skew- symmetric 

(e-a,b) = -(a,e-b) (15) 

with respect to the scalar product 

/ + 00 
f (a;) ■ g(x) dx , (16) 

-co 

and if it satisfies a "Jacobi-like" identity. The precise form (A.l) of that identity is given in 
the Appendix, since it is not of immediate importance here. A vector function 7^(-u,f) 
is a gradient function if its Frechet derivative is symmetric with respect to the scalar 
product (16): 

(7^[a],b) = (a,7^[b]) . 



Recall that the Frechet derivative of a vector function 7^(zx,t;) 






(17) 



m a 



direction 



V 



is given by 



7h 





e=0 



7i(m + e^, v) + 7i(m, v + er]) 
l2{u + £^, v) + 72(m, V + er]) 

The Hamiltonian, or Hamiltonian functional, H = J^^ p{u, v)dx giving rise to the gradi- 
ent function ^^{u.,v) through the identity (Berger 1977) 



J —OD 



(19) 



can be recovered from this gradient in the following way: 



H 



u 



{-f^{\u,\v),\ \)-C 

V 



d\ , 



(20) 



where the constant C is chosen such that the integral giving H is bounded. Such a 
Hamiltonian formulation admits a Poisson bracket, 



{A 5} = (7a, e- 7b) 



defining the time evolution of a functional A oi u and v. 



(21) 



dA _ dA a' I ^ 
dt dt \ V 



using formulas (19) and (14). 



dA , 



u 

V 



dA 
'dt 



+ {1a.Q-1h) 



dA 
'dt 



+ {A^}, 



(22) 



Since O is skew-symmetric, and due to identity (A.l), the bracket (21) possesses all 
the characteristics of a standard Poisson bracket, except for the "Leibniz-like" expul- 
sion property which cannot be properly defined in the case of functionals. Clearly, if 
autonomous, the Hamiltonian H itself is a conserved quantity for the evolution equation 
(14). Since (6) can be written as a conservation law, an apparent Hamiltonian formulation 
is the following: 



u 

V 



02-72> 



(23) 



dr 

a 



©2 = - I „ „ I , 72 

where dx denotes the partial derivative with respect to x 



{u^ + v'^)v + Ux 



O2 is easily seen to be a Hamiltonian operator since it is a constant {i. e. not depending 
on u or v), skew-symmetric operator (see Appendix). 

It is straightforward to verify that 72 satisfies (17) and thus leads to a Hamiltonian 
functional H2 given by formula (20): 



Ho 



dX 

j-00 



X\u^ + v^y + X{uxV - uvx + (3u^) - C2 



dx 



(24) 



{u^ + V^Y UxV — UVx P 



+ 



+ '-u^ - C2 



dx . 



Comparison with formula (10) shows that it is the conserved density p2 which give rise 
to this Hamiltonian. One may wonder if the conserved density pi in (9) can be hnked to 
a Hamihonian functional as well. Let us define Hi by 

n+oo r , 



Hi 



Pi ^ 



dx 



u^ + f ^ 



2 
Using (19), the gradient 7^^ of this functional is found to be: 



-Ci 



dx . (25) 



7i = r 1 . (26) 

V 



The challenge is to find a Hamiltonian operator ©i such that (6) can be recast into the 
form (14) using 7^^. The MVDNLS equation would then have a second Hamiltonian 
formulation and thus possess a bi-Hamiltonian structure. 



5 Bi— Hamiltonian structure and recursion operator 



If an evolution equation admits a bi-Hamiltonian formulation (Magri 1978) 



u 



V 

■ t 



ei -71 = 62 -72, (27) 



and if it is possible to show that Oi and ©2 are compatible Hamiltonian operators {i.e. 
©1 + ©2 is again a Hamiltonian operator), and if one of the operators, say ©2 is invertible, 
then the operator 

R = Q^. ©-1 (28) 

is a hereditary recursion operator (Fuchssteiner 1979, Fuchssteiner and Fokas 1981) for 
that evolution equation. 

The formal adjoint of R with respect to the scalar product (16) 

i?t = 0-1 . 0^ (29) 

then maps gradients of conserved quantities into gradients, provided the operator R is 
injective (Fuchssteiner and Fokas 1981, and also Appendix), thus defining an infinite 
sequence of conserved quantities, all in involution with respect to the Poisson bracket 
(21). 

In our case we have already found one Hamiltonian structure with an invertible Hamil- 
tonian operator ©2: 



where 9^ ^ denotes the inverse of the d^ operator, such that d^d^ ^ = d^^d^ = 1- 

A first step towards finding a second Hamiltonian structure is the construction of 
a skew-symmetric operator ©i which casts (6) into Hamiltonian form with the gradient 
7i in (26). The most general parametrization (involving a single d~^ operator) which 
satisfies the above constraints is 

e, ^ - ( '; '; ] , (31) 

y f2i fc'22 J 

with 

6*11 = (idx + civd^v + csudxU + (6 - 4c3)u^d~'^Ux , 

a a2 I / 3 Ci + C2 C5 o , / 1 Ci — C2 C5 

a «2 , /3 C1 + C2 C5 1 C1-C2 C5 ^ ^ 

^21 = ^. + (2 2 X^ "^2 ^~ T^^"^ 

+ (1 - ci - -^)uv^ + C5t;^a~^M^ , 
6*22 = C2udxU + Civd^v + (6 - Ac^jv^d'^v^ . 

At this point one could of course try to find values of the parameters for which this 
operator is actually a Hamiltonian operator. However, this proves to be a formidable task 
(see Appendix). A better approach is to concentrate on the action of R'^ on the gradient 
functions obtained so far. Because of the bi-Hamiltonian structure (27) of the equations, 
B) maps 7i into 72: 

72 = e2-i-0i-7i = i?^-7i. (33) 

Suppose R'^ actually is the formal adjoint of a recursion operator for the MVDNLS equa- 
tion, then it will map 72 into yet another gradient function. If the conserved quantity 
corresponding to this gradient has to be polynomial in u, v and their derivatives, then the 
values of the parameters in (32) have to be such that ©i ■ 72 is a total x-derivative: 

73 = e2-^-01-72 = -rf g\]-®^-^2- (34) 

This requirement uniquely determines all five parameters in (32), namely Ci = C2 = 0, C3 = 
C4 = — C5 = 2, thus giving us a single candidate for a second Hamiltonian operator: 

9 ^ _ ( /^"^x + 2ud^u - 2ua;d-^u^ -dl + 2vd^j,u - 2u.j,d-^v^ . 
\ dl + 2udxV - 2vxd~'^Ux 2vd^v - 2v^d~'^v^ 

together with the gradient 

^' ~ ' \ -Pv^ - V,, + Iv' + lu^v + 3u^v^ + 3u,v^ + 3u^u, ' ' ^ ' 

10 



Redefining 73 as 73 — 13^ 21 i^ i^ straightforward to sliow tliat it satisfies (17) and corre- 
sponds to the Hamiltonian functional 



H. 



{U^ + V'^f 3 3 1 , . , P f A 4n ^ 

Vu^v -u v^- -[uu^^ + VV^^) + j[u -V ) -C3 



dx . (37) 



The density function associated with this functional is, up to partial integration, the 
conserved density ps given in (11). Hence, it follows that i?^ maps 72 i'^to the gradient of 
a conserved quantity (i.e., H3 + PH2), suggesting that R indeed is a hereditary recursion 
operator for the MVDNLS equation. 

In the Appendix it is proven that ©i is a Hamiltonian operator and that it is com- 
patible with 02. Thus, we have shown that 

^ ^ / /3 + 2^2 + 2u,d-^u -d, + 2uv + 2u,d-^v \ 
~ \d^ + 2uv + 2v^d-^u 2v^ + 2v^d-^v ) 

which follows from (28), (30) and (35), is a hereditary recursion operator for the MVDNLS 
equation: 

The recursion operator therefore defines a hierarchy of integrable evolution equations 






= R--\ I , (40) 

all sharing an infinite sequence of conserved quantities, the gradients of which are gener- 
ated by the formal adjoint of R: 

^t^( P + 2ud-'ud. -d + 2ud-'vd, \ ^^^^ 

\ dx + 2vd~^udx 2vd~^vdx j 



6 The DNLS as special case 



For /3 = we recover from the MVDNLS the usual DNLS itself. Hence R\p=Q or 

^^( 2u' + 2uxd-^u -dx + 2uv + 2uxd-^v \ 
~ \dx + 2uv + 2vxd-^u 2v^ + 2vxd-^v ) 

is the recursion operator for DNLS, since it can be seen that the /? term in Oi does not alter 
the fact that it is a Hamiltonian operator (see Appendix). The Hamiltonian formulation 
(23) reduces to the one given by Kaup and Newell (1978) when /J = 0. However, to the 
best of our knowledge, the second Hamiltonian formulation is a new result, not found in 
the literature, and neither is the recursion operator (42). 

11 



7 Conclusions 

Slightly oblique propagation of Alfven waves in strongly magnetized plasmas is described 
by a nonlinear vector evolution equation which differs from the vector form of the DNLS 
by an extra linear term, if one retains both the deviation from charge neutrality between 
the different plasma species and the influence of the displacement current in Ampere's 
law. The resulting modified vector derivative nonlinear Schrodinger equation cannot be 
transformed into the DNLS itself, and this has implications for the discussion of integra- 
bility and the possibility of finding conserved densities and solitary wave solutions. Of 
course, the reverse is true: the MVDNLS includes the DNLS as a special case. 

While in a previous paper the applicability of the prolongation method, adapted 
to a vector nonlinear equation, and the existence of some invariants were indicative of 
complete integrability, in the present paper we have shown that the MVDNLS indeed 
possesses a bi-Hamiltonian structure, and hence, through the resulting recursion operator, 
an infinite sequence of conserved densities. We were guided in this by the explicit symbolic 
computation of the first seven conserved densities. 

Surprisingly enough, there seems to be no proof in the literature that the DNLS itself 
has a bi-Hamiltonian structure, although the existence of an infinite sequence of conserved 
densities was known. As the DNLS emerges as a special case of the MVDNLS studied 
here, we now also have the explicit form of the recursion operator for the DNLS, besides 
the proof of its bi-Hamiltonian structure. 
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Appendix 

A skew-symmetric operator is a Hamiltonian operator if and only if it satisfies the 
identity (Fuchssteiner and Fokas 1981) 



(a, 0'[0 ■ c] ■ b) + (b, 0'[0 ■ a] ■ c) + (c, 0'[0 ■ b] ■ a) = , 



(A.l) 



where the scalar product was defined in (16). Since this identity involves taking the 
Frechet derivative of in a certain direction, it is trivially satisfied if does not depend 
upon u or v, i.e. if it is a "constant" operator. Hence, every constant skew-symmetric 
operator is a Hamiltonian operator. 

Verifying that 0i satisfies (A.l) is a very cumbersome task. A method which makes 
such a verification a lot more tractable relies on defining a functional multivector (see 
Olver 1986) 



^4/{(c.)Ae-(^)> 



dx . 



(A.2) 



Using this formalism, it can be shown (Olver 1986) that the identity (A.l) can be recast 
into the form 

V 
assuming, by definition, that C, and rj are independent of u, v or their derivatives. 



Z 







= 0, 



(A.3) 



In our case we define the multivector 



+ AuxV^ Arj — 4:Ux^ A d~^Vxi] + 2v'^rj Arfx — 2f^7] A d'^Vxt]} dx . 
It is a straightforward, albeit tedious, calculation to show that 



(A.4) 



Z, 



Qi 



7] 



0, 



(A.5) 



thus proving that 0i is a Hamiltonian operator. 



In the same manner it can be shown that 



Z[ 



0. 



V 



0, 



(A.6) 



and hence that 0i + 02 is a Hamiltonian operator as well, i.e. thus 0i and 02 are 
compatible Hamiltonian operators. Note that since the only term in Zi depending on (3 
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is independent of u or v, its presence (or absence) does not alter the properties (A. 5) or 
(A. 6). Consequently, the above considerations also hold in the special case of the standard 
DNLS equation. 

Concerning the injective nature of the recursion operator R, which is needed to guar- 
antee that W maps gradients to gradients (see Proposition 2 in Fuchssteiner and Fokas 
1981), it suffices to consider the action of R on the space Sp of (2-component) vector 
functions with polynomial components in u, v and their derivatives, which are all C°° 
functions fast decreasing at infinity to symmetric but in the case of u non-zero boundary 
conditions. It is easily seen that the kernel of R in this space consists of the elements 
which are of zero polynomial degree. Since R either raises the polynomial degree of the 
different terms by two or leaves it invariant, R can be made injective by restricting the 
domain to the space Sp^ C Sp without zero degree elements. 
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